Excitation spectra of the two dimensional Kondo insulator 
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We present a Quantum Monte Carlo (QMC) study of the temperature dependent dynamics of the 
two-dimensional (2D) Kondo insulator. Working at the so-called symmetrical point allows to per- 
form minus-sign free QMC simulations. Study of the temperature dependence of the single-particle 
Green's function and the dynamical spin correlation function provides evidence for two characteristic 
temperatures, which we associate with the Kondo and coherence temperature, whereby the system 
shows a metal-insulator transition at the coherence temperature. The data shows evidence for two 
distinct types of spin excitations and we show that despite strong antiferromagnetic ordering at low 
temperature the system cannot be described by spin-density-wave (SDW) theory. 
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The periodic Anderson model (PAM) or its strong 
coupling version, the Kondo lattice, may be viewed as 
the appropriate model for describing intensively investi- 
gated classes of materials as the heavy electron metals 
and the Kondo insulators [^IJ]. While the impu- 
rity case is well-understood |^-^ and even amenable to 
exact solutions [|lO|,0, little is known about the lattice 
model. Recently a considerable amount of numerical re- 
sults has been collected, albeit mainly for the the one- 
dimensional (ID) case ]T2| , p^ . The main reason is that 
the frequently used density matrix renormalization group 
(DMRG) method works best for ID systems. Only very 
recently finite-temperature exact diagonalization results 
for the two-dimensional (2D) case | |T^ b ecame available. 

on the other 



The Quantum Monte-Carlo method |lj-|ri_ 
hand, can in principle treat systems of arbitrary dimen- 
sion. Here the limitation is in the notorious minus-sign 
problem, which usually precludes the study of truly low 
temperatures. However, by restricting oneself to the so- 
called symmetric point the minus sign problem can be 
circumvented and quite low temperatures be reached in 
numerical simulations |18 1^. In the present manuscript 
we want to present data for the dynamics of the 2- 
dimensional (2D) PAM with Hamiltonian 



H 



H.c. 



(1) 



Here cj ^ {f} creates a conduction electron (/-electron) 

in cell «, rii^cT—fl afi a ^^^^1 denotes nearest-neighbor 
lattice sites. The so-called symmetric case corresponds 
to the special choice U — 2ef. At 'half-filling', i.e. elec- 
tron density /unit cell n — 2, the Hamiltonian acquires 
particle-hole symmetry, whence there is no more minus 
sign problem. 



In a preceding work |19| we studied the symmetric ID 
case with interaction value U = 8.0 1 and c-f hybridiza- 
tion V = 1.0 1. For the lowest accessible temperature 



(T — 0.03 less than 1% of the conduction-electron 
bandwidth) the systems exhibits insulating behavior with 
a gap in the single-particle spectrum. The c- and f- 
electrons seem to form a coherent all-electron fluid with 
composite c-f character of the low-lying one- and two- 
particle excitations, thus turning the system to a 'nom- 
inal' band insulator with a 'Fermi-Surface' which covers 
the entire Brillouin zone. Above the lower crossover tem- 
perature Tcoh (which we identified with the so-called co- 
herence temperature in heavy Fermion systems) the c- 
and f-like features in the correlation functions are de- 
coupled, indicating that the f-electrons 'drop out' of the 
Fermi-surface. The Fermi surface shrinks to one half of 
the Brillouin zone whence the system becomes a metal 
and the single-particle gap closes. At Tcoh the spin ex- 
citations of the f-electrons become localized (visible as a 
dispersionless branch in the f-electron spin response) and 
the spin gap closes due to c-like spin excitations. Kondo 
resonance-like sidebands with f-character are seen in the 
single-particle spectral function presumably as a sign of 
the formation of loosely bound singlets between c- and f- 
clectrons. At the higher crossover temperature Tr- these 
dispersionless f-like Kondo resonance-like sidebands in 
the single-particle spectrum disappear, whence we iden- 
tify it with the Kondo-temperature of the system. Above 
this temperature the single-particle spectral function of 
the c-electrons shows a very conventional tight-binding 
dispersion — 2tcos(A;), whereas the f-electrons show nor- 
mal upper and lower Hubbard bands, i.e. the f-electrons 
do not participate in the low-energy physics. While this 
is hard to establish numerically, it follows from general 
theorems for ID systems that there is no long range order 
at any temperature. 

This does not appear to be the case in the 2D case stud- 
ied previously with a smaller interaction U — 4.0 i by 
Vekic et al. iQ. There the ground state of the sys- 
tem is an insulator with long-range AF order, a finite 
charge gap and gapless spin excitations for small values 
of V^/Ut (i.e. a Mott insulator). As V^/Ut increases. 
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the long-range order is destroyed and spin-liquid behav- 
ior is found characterized by both a spin gap Ag and a 
charge gap Ac with Ac > A5 (i.e. a Kondo insulator). 
Further increasing the hybridization V then even leads to 
a band-insulating state with equal values of the spin and 
charge gap. The authors also found a Kondo resonance- 
like peak in the angle-integrated spectral density 
for moderate temperatures and a gap for temperatures 
below. 

Similar results concerning the ordered nature of the 
ground state in 2D were obtained in a recent zero tem- 
perature QMC study of the strong coupling version of the 
model by Assaad |2Cf| . There, a quantum phase transition 
was found to take place from an antiferromagnetically or- 
dered phase for J/i < 1.5, to a spin liquid phase for larger 
J/t. Here J denotes the Kondo exchange which would 
come out as J = 8V'^/U in the strong coupling pertur- 
bation approach due to Schrieffer and Wolf [|T]. One 
would thus expect that in the Hamiltonian (|l|) the tran- 
sition occurs at V'^/Ut w 0.18 in reasonable agreement 
with the results of Vekic et al. for U=4. 



filling, i.e. with two electrons/unit cell, again at the sym- 
metric point ej = U/2 (i.e. J/t=l). Due to the absence 
of the minus-sign problem we could reach temperatures 
as low as T = 0.05 1, corresponding to w 0.6% of the 
conduction-electron (c-electron) bandwidth. 
Based on the previous works ||l^,^ we expect that due 
to our larger value for the interaction U = 8.0 1 with a 
smaller ratio of V'^/Ut, an insulating ground state with 
long-range AF order to be stable. This is confirmed by 
our numerical results for the spin and charge susceptibil- 
ities Xs Xc = c, / denotes the type of electron 
probed). These are defined as 



Xc 



dM 
dN 



(n"(T)n"(0))dT, 



(2) 



where 5'^ = i(n^ -nj), = (n^ -Hnp. Figures and | 
show these susceptibilities for the ID and 2D lattice. In 
both cases the charge susceptibility for the conduction 
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FIG. 1. Susceptibilities of the ID Kondo lattice with 
U = 8.0 1 and V = 1.0 1 for different temperatures: (a) Charge 
susceptibility Xcc, (b) spin-susceptibility Xaz and (c) inverse 
spin-susceptibility l/xsz- 

In this work we study the 2D Kondo lattice with 6x6 
unit cells using standard QMC techniques for the same 
interaction U = 8.0 1 as in our previous work for the ID 
case Q . Again we restrict ourselves to the case of half 
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FIG. 2. Susceptibilities of the 2D Kondo lattice with 
U = 8.0 1 and V = 1.0 1 for different temperatures: (a) Charge 
susceptibility Xcc, (b) spin-susceptibility Xaz and (c) inverse 
spin-susceptibility 1/Xsz- 

electrons are considerably larger than those for the /- 
electrons. This simply reflects the fact that the /- 
electron subsystem is practically half-filled and the strong 
Coulomb repulsion renders the /-electron system incom- 
pressible. At high temperatures Xc the ID system 
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is more or less temperature independent, whereas it in- 
creases with decreasing temperature for the 2D system. 
This can be understood by assuming that in both cases 
/ and c-electrons are decoupled at high temperature, as 
is suggested by our previous study for the ID model [ p^ . 
The different temperature dependence then can be traced 
back to the different free-electron density of states (DOS) 
in ID and 2D. Namely for free electrons with /i(T) — 
one has Xc ~ where fc^ is the Boltzmann constant 
and p is the DOS averaged over a window of width fc^T 
around the chemical potential jj, = 0. This will be more 
or less constant for ID, but strongly increasing with de- 
creasing T in 2D, where one has a van-Hove singularity in 
the band center. At temperatures below O.lOt, however, 
Xc drops sharply in both ID and 2D. This is the behavior 
expected for a band insulator and indicates the opening 
of a gap in the fully interacting DOS, because shifting /i 
does not change the particle number any more. In both 
ID and 2D we thus have evidence for the opening of a 
gap in the single-particle spectrum. 



A related quantity of interest is the static magnetic struc- 
ture factor 
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FIG. 3. Antiferromagnetic structure factor Saf as a func- 
tion of temperature for /-electrons (full line) and c-electrons 
(dashed line), calculated for the 2D 6x6 cluster (open 
squares) and a ID 16-cell chain (dark squares). 

Turning to the spin susceptibility we note that there the 
dominant contribution comes from the /-electrons. In 
both ID and 2D xl can be fitted roughly to a Curie- 
law at high temperatures. Deviations occur for smaller 
T < 0.20 i. It is also in this temperature range that the 
difference between the ID and 2D system becomes ap- 
parent: the downturn of the spin-susceptibility XsziT) at 
the coherence temperature Tcoh for the ID case signals 
a ground state with a spin-gap and thus without long- 
range AF order, whereas the steady increase of the 2D 
spin-susceptibility Xsz{T) suggests a ground state with- 
out a spin-gap and thus with long-range AF order in 
agreement with the data of Vekic et al. |l8j and Assaad 
pO| . We also note the overall agreement of the suscepti- 
bilities with the exact diagonalization data by Haule et 



S{q) = {Sl^S:^) 



(3) 



Figure |^ shows the antiferromagnetic structure factor 
Saf, i-e. S{q) for g = tt in ID and for q = (tt, tt) in 
2D. At high temperatures, Saf is T-independent and for 
the /-electrons is close to the value of 0.75 expected for 
a totally uncorrelated spin-I/2 system. Saf for the c- 
electrons is equally T-independent, but takes a smaller 
value presumably due to the reduction of the c-moment 
by the charge fluctuations of a free-electron gas. At low 
temperature Saf for the /-electrons in 2D increases quite 
dramatically, indicating the tendency towards long-range 
order. The c-electron structure factor on the other hand 
shows only a weak enhancement at low T: the order- 
ing seems almost exclusively restricted to the /-electrons. 
The Saf in ID also increases, but seems to saturate at 
lower temperature. In our preceding study in ID we 
found a spin correlation length of C w 4.6 at the lowest 
accessible temperature T = 0.03 i in a ring of 16 unit 
cells. The fact that ^ was already quite significantly 
smaller than the cluster size indicated that long range 
order does not develop in this case. A similar analysis in 
2D is difficult, because the clusters we can study have a 
much smaller linear extent than in ID, but it is visible 
that the two systems behave quite different in their mag- 
netic properties at low temperatures. As we will see in a 
moment, this difference at low temperatures has practi- 
cally no bearing for the temperature dependent dynamics 
of the models, which are practically indistinguishable in 
ID and 2D. In particular we will see that the rather anti- 
ferromagnetic correlations in the 2D system do not affect 
the single-particle spectra in any noticeable way. 
We start in Figure ^ with a discussion of the tempera- 
ture development of the angle-integrated spectral density 
of states D{ijj) as a function of temperature (see Ref. |19| ] 
for a definition of the various correlation functions). At 
the highest temperature we studied, T = 1. 00 1, the f- and 
c-electrons are completely decoupled and for both ID and 
2D system we observe nearly identical f-like upper and 
lower Hubbard bands at « ±?7/2, whereas the c-electrons 
show the density of states expected for free electrons in 
the respective dimension. This decoupling of the c- and 
the f-electron physics is also visible in the momentum re- 
solved single-particle spectral function A{k, cj) plotted in 
Figure || (now and for the rest of this work only for the 
2D case; the similarity between ID and 2D is nevertheless 
striking and we refer the readers to the plots in our pre- 
ceding publication |I9|). At temperature T = I.OOt the 
c-electrons show a very conventional free tight-binding 
dispersion, whereas the f-electrons obviously do not par- 
ticipate in the low-energy physics at all. The momentum 
resolved spin-correlation function S{k,uj) (see Figure^) 
shows a free-electron-like particle-hole continuum in the 
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c-electron case and a practically dispersionless branch in 
the f-like spectrum, identifying the magnetic f excitations 
as practically immobile. 



0.4 
0.3 
0.2 
0.1 


0.3 
0.2 
0.1 


0.3 
0.2 
0.1 


0.3 
0.2 
0.1 


0.4 
0.3 
0.2 
0.1 





- T^l.OOt 


f-electron 

c-electron 

Id f-elec. — 

Id c-elec. 


r = 0.3 






T = 0.20r 




r = 0.10? 




T = 0.05r ), 

_ I J 







(O-fj/t 



8 



FIG. 4. Angle-integrated spectral density D{u!) of the 2D 
Kondo lattice with U — 8.0 1 and V — 1.0 1 for different tem- 
peratures. The ID Kondo lattice density of states D{uj) is 
shown at temperature T — l.OOt for comparison. 

Lowering the temperature to T = 0.33 1 the density of 
states D{uj) shows f-like spectral weight around /i which 
forms a single peak right at fi as the temperature is 
lowered further to T = 0.20 1. We interpret this as a 
sign of the formation of loosely bound singlets between 
f- and c-electrons, i.e. the /-electrons now start to par- 
ticipate in the low-energy physics. At this temperature 
the c-electrons also start to deviate from the free tight- 
binding dispersion with a broadened and slightly split 
peak right at /i for momentum k — (tt, 0) and with repli- 
cas of the tiny f-like "foot" , best visible at the neighbor- 
ing momenta k = (37r/4, 0) and k = (7r,7r/4), but also 
for k = (7r/2,7r/2). At T = 0.20t the dynamical spin- 
correlation function S{k,uj) of the f-electrons shows still 
a dispersionless branch, i.e. the magnetic f excitations 
are still immobile, but now with a much more narrow 
width. In addition broad c-like low-energy 'humps' are 
being split off from the c-like electron-hole continuum, 
best visible along (0,0) (7i',7r). 



Proceeding to the temperature T = O.lOt the f-like den- 
sity of states in Figure |4| shows the formation of side 
bands at « ±li with the c-like peak right at /i starting 
to split as precursor of the formation of the single-particle 
gap. 
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FIG. 5. Momentum resolved single-particle spectral func- 
tion A{k, uj) with same parameters as in Figure Wifor different 
temperatures. 

At the lowest temperature T — 0.05 1 the c-electron den- 
sity at higher energies is still consistent with a standard 
2D tight-binding density of states (the maxima at « ±2 1 
are a finite-size effect - we have checked that they appear 
also when a free tight binding band is simulated in a fi- 
nite cluster). At low energies, however, D{uj) shows a 
clear gap around fi. This demonstrates the insulating 
nature of the ground state, in agreement with the be- 
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havior of the charge susceptibihty Xcc- The f-electrons 
show the high-intensity upper and lower Hubbard bands 
at w ±C^/2, but now also very sharp low-energy peaks 
at the gap edges and in addition the two side bands at 
K, ±1 1. These sidebands are best visible in the spectral 
function A(fc, uj) of the f-electrons shown in Figure |^ and 
are accompanied by a further change in the spectral func- 
tion of the c-electrons: The cos(fca;) -|-cos(fcj,) dispersion 
of the c-electrons now shows a gap at fc^ — (tt, 0). 
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FIG. 6. Momentum resolved dynamical spin-correlation 
function S{k,uj) with same parameters as in Figure ^for dif- 
ferent temperatures. 

In the 6x6 cluster this is the only momentum on the 
Fermi surface for noninteracting conduction electrons, 
but one may expect that in a larger system the gap is 
uniform in fc-space. The states with the lowest exci- 



tation energies now are practically dispersionless bands 
with small, /-like weight. Obviously the /-electrons now 
fully participate in the single-electron states close to fi. 
Assuming the validity of the Luttinger theorem with the 
number of electrons being given by c- and /-electrons 
together would also give an obvious explanation for the 
insulating nature of the ground state - at electron density 
n = 2 the system should be a band insulator. With the 
exception of the /-like side bands the single particle spec- 
trum agrees very well with theoretical predictions based 
on an 'expansion around an singlet vacuum' p^ . 
The dynamical spin-correlation function S{k, uj) of the 
f-electrons at T = 0.05 t shows an intense branch of 
low-energy excitations with a weak dispersion. Its spec- 
tral weight has a sharp maximum at fc = (7r,7r), consis- 
tent with the strong AF correlations for the /-electrons. 
In the c-electrons' spin correlation function there ap- 
pears a practically dispersionless low energy excitation at 
u! « 0.50t. In some cases one can see that the /-electron 
spectrum also shows a weak hump at the position of this 
dispersionless excitation, indicating that this excitation 
has a mixed f-c character. The sharp /-like low energy 
mode also has some admixture of c-weight at (tt, tt) - the 
spin excitations at low energies thus have both / and c 
character, which shows again that the two types of elec- 
tron have merged to form a common all-electron fluid. 
As a surprising result the temperature development of 
the dynamical correlation functions resembles in con- 
siderable detail that seen previously in ID - even 
the 'crossover-temperatures' in the spectral function are 
practically the same in ID and 2D. In fact, with the sole 
exception of the different form of the c-electron DOS at 
higher energies, the angle integrated spectra in Figure ^ 
are practically indistinguishable from their ID counter- 
parts in Figure 1 of Ref. ||l9| , and this holds true for each 
individual temperature. Assuming that this is not coin- 
cidence for the special set of parameters we are using, 
the characteristic temperatures for the ID and 2D mod- 
els with identical parameters thus are at least very close 
to one another. This would be hard to understand if we 
were to assume that the characteristic temperatures of 
the model depend sensibly e.g. on the c-electron density 
of states at the Fermi energy, pp - this would be singular 
(or at least significantly larger on a discrete lattice) for 
the 2D case. In the case of the Kondo-impurity, where 
the Kondo temperature is Tk oc exp{—j^), we would 
thus expect a very different temperature evolution. To 
be more quantitative, we estimated the impurity model 
Kondo temperatures for finite systems. To better take 
into account the particle-hole symmetry we assume that 
the transformation to the strong coupling model has been 
performed, i.e. we use 
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H 



(4 



H.c.) + J^S ■ 



,Ck' 



k,k' 



(4) 



where t denotes the vector of Pauh matrices, and S 
the spin operator for the /-hke impurity. The value of 
J — 8V'^ /U thereby takes into account the fact that the 
exchange processes can occur both via an empty or a dou- 
bly occupied /-level. We then use the variational trial 
state by Yoshimori which yields a self-consistency 
equation for the binding energy Es of the Kondo singlet: 



1 



3J 



1 



4N ^ Ek - Es 

K^unocc. 



= 0. 



(5) 



To solve the equation, we take the true discrete fc-meshes 
of the 16-site chain and the 6x6 lattice, which presum- 
ably takes finite-size effects into account to some degree. 
We then find for our parameters the singlet binding en- 
ergies Es = -0.0624 t in ID and -0.1219 t in 2D. The 
Kondo temperatures in the impurity model thus would 
differ by roughly a factor of 2, but in the lattice model 
the characteristic temperatures are more or less indis- 
tinguishable. Assuming that the (near)-identity between 
the characteristic temperatures in ID and 2D is not just 
a coincidence for the specific parameter set we are study- 
ing, this suggests that the characteristic temperatures for 
the lattice system have little or no relationship with those 
for the impurity model. 
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FIG. 7. Momentum distribution for /-electrons (full 
squares) and c-electrons (light squares) at T = 0.33t (full 
lines) and at T = O.OSt (dashed lines) along high symmetry 
lines of the Brillouin zone for the 6x6 cluster. 

A further surprising feature of the results is the absence of 
any sign of antiferromagnetism in the single-particle spec- 
tra: the fe-resolved spectra for the c-electrons do not show 
any indication of antiferromagnetic umklap bands (for 
the /-electrons it is impossible to make such a statement 
because the /-like bands are all more or less flat). Figure 



shows the momentum distribution ~ X^o-i'^I a^k,<j} 
(a = c, /) at high temperature T = 0.33i, where no 
noticeable enhancement of Saf could be seen in Figure 
^ and at the lowest temperature T = 0.05i where Saf 
shows strong signatures of antiferromagnetism. The most 
prominent feature is the 'Pseudo Fermi surface' for the 
c-electrons, i.e. a sharp drop of which occurs at the 
Fermi surface of the unhybridized conduction electrons. 
This is familiar from numerical studies of the ID model 
| p2| , p^ and can be explained theoretically by 'expansion 
around the singlet vacuum' [ p2| . Surprisingly, is al- 
most indistinguishable, the only change at low tempera- 
ture being a 'sharpening' of the distribution near (tt, 0). 
Similarly, develops some structure at low tempera- 
ture and is completely flat at high temperature (the lat- 
ter again shows the decoupling of c and /-electron at 
high temperatures). We note that both changes are in 
fact opposite to what one would expect on the basis of 
a conventional spin-density-wave (SDW) picture: there, 
the static SDW would provide an additional potential 
V{Ri) = VAFe''^ '' which would tend to mix the single 
particle states a\. ^ and Q;^_|_q ^. Any difference between 
in the inner part and the outer part of the antifer- 
romagnetic Brillouin zone should thus be reduced by the 
antiferromagnetic ordering. By contrast, the actual data 
show that the structures in actually sharpen up in the 
ordered state. 

To be more quantitative, let us briefly discuss inhowmuch 
an SDW- like mean- field theory might explain our results. 
We approximate the interaction term for the /-electrons: 



(6) 



Here Uf is the average density of /-electrons/site, 
and particle-hole symmetry at half-filling implies 
nf=l. The parameter m is the staggered magneti- 
zation of /-electrons. Introducing the vector C — 
ick,cr,Ck+Q,crfk,a, fk+Q,a), thc Hamiltouiau then takes 
the form 



J2 Y.cLHk.^c^ 



(7) 



keAFBZ a 



with the matrix 



efc, 0, 

0, -tk 

-V, 0, 

0, -F, 



-F, 
,0, -V 

2 

-(TmU fx 



(8) 



For our parameter values, self-consistent calculation of 
the staggered magnetization m yields the value m « 0.95. 
The band structure obtained in this way is shown in Fig- 
ure ^. It can be understood by considering the limit 
V/U ^1 where we expect first of all dispersionless /-like 
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bands at ±17/2 (obtained by diagonalizing the 2x2 ma- 
trix in the lower right corner of (H)). EUminating these 
bands from the Hamiltonian by canonical perturbation 
theory generates a coupling matrix element between Cfc^o- 
and Ck+Q^a- which is equal to —2V'^/U. We thus expect 
for the c-electrons a band structure which resembles the 
SDW-approximation for the single-band Hubbard model, 
with the relatively small SDW gap parameter 2V'^ jU . 
The dispersionless /-like bands close to /i, which can be 
seen rather clearly in the numerical spectra in Figure ^, 
are not at all reproduced by this approach, whence a sim- 
ple SDW-like mean-field theory is clearly inadequate to 
explain the band structure generated by the QMC simu- 
lation. 
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FIG. 8. Self-consistent SDW band structure for the same 
parameter values as Figure ^ 

The above results allow some conclusions concerning the 
spin excitations responsible for the ordering if we consider 
the possible states of a single cell. If we restrict ourselves 
to low-energy states, the most likely states are those with 
precisely one /-electron/unit cell. These would be 

Kl,a)-/2J0), 
K3,a)=4^c[^4j0), 

K-2,0) = ^(<,4-ct,/t^)|0), 

N,2,l) = -^(4,/^^ -1-4,4^)10). (9) 

These states have a spin of 1/2 cr) and |«,3, cr)), 

(I?, 2,0)) or 1 (|i,2,l)). Above, we have defined the 
5"^ = component of the triplet, but there are also states 
with Sz — ±1. 

Then, one can envisage two very different types of spin 
excitations: 

a) One could flip the spin of one of the charged cells, e.g. 
convert |1, ct) — > |3, ct). In these states the spin of the cell 
is carried exclusively by the /-electron. This means that 



this type of spin excitation can be created only by acting 
with the /-spin operator. The c-electrons' spin operator 
obviously cannot 'touch' these types of single-cell states, 
b) Alternatively, one could convert the singlet into one 
of the three components of the triplet. In this case, the 
total spin of the cell is carried by both, the c and the 
/-electron. This type of excitation therefore can be gen- 
erated by both, the /- and the c-spin operator. 
Since the two types of spin excitations have identical 
quantum numbers, one might expect that the 'true' spin 
excitations of the system are a mixture of the two. How- 
ever, our data suggest that the mixing between the two 
types of excitation is indeed quite weak: for example the 
fact that it is predominantly the /-like structure factor 
Saf which shows ordering in Figure ^ indicates that the 
AF ordering is predominantly due to ordering of the /- 
spins in singly and three-fold occupied cells. Let us for 
example define 

K a) = cos(e)K 2, 0) + ^^(N, 1, ^) + K 3, a)), (10) 

and (introducing the two sublattices A and B) 

m=v2NX{\i,])X{\iA), (11) 

where 'P2N projects onto states with precisely 2N elec- 
trons and z-spin Sz — 0. The state j^*) then has Neel or- 
der with an ordered moment cx sin^(0) in the /-system, 
but no order whatsoever in the c-electron system. We 
do not claim that the state l^t^) has much to do with 
the ground state of the lattice model - it demonstrates, 
however, that by using singly and three-fold occupied 
cells it is indeed possibly to construct states where only 
the /-electrons do order. This suggests that the order- 
ing is driven by the charged single-cell states |i, 1, cr) and 
\i,3,a) (which should be modeled as effective Fermions 
p2| ) and not the singlet-triplet excitation. This assump- 
tion also explains immediately ,why the c-electrons mo- 
mentum distribution (see Figure ^ and spectral density 
(see Figure |^) is so remarkably unaffected by the AF or- 
der at low temperature. 

Rather clear evidence for these two types of spin exci- 
tations is also provided by the dynamical spin correla- 
tion function at low temperature. There, we have seen 
at the lowest excitation energies an almost purely /-like 
spin excitation with some relatively small admixture of c- 
weight near (tt, it). This would suggest that this mode has 
the character of a particle-hole excitation carried by the 
charged spin- 1/2 cells, |i,l,cr) and |i,3, u). At a some- 
what higher excitation energy, the data showed a practi- 
cally dispersionless mode with strong c-character but also 
some /-character - this might correspond to the singlet- 
triplet excitation. The weak admixture of c-weight at 
(tt, tt) in the low energy spin excitation then is a measure 
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for the mixing between the two types of spin excitation. 
In summary, we have presented a QMC study of the tem- 
perature dependent dynamics in the 2D Kondo lattice 
model at half- filling. Working at the symmetric point, 
U = 2e/, we could avoid the minus-sign problem and 
study the evolution with temperature down to extremely 
low temperatures. As was the case in ID [Q, we could 
identify two characteristic temperatures. The lower one 
of these, which we identify with the experimental coher- 
ence temperature, is associated with the metal-insulator 
transition in the half-filled case, n = 2. Below Tcoh the /- 
electrons merge with the conduction electrons to form a 
coherent all-electron fluid. The /-electrons consequently 
participate in the Fermi surface volume, whence the sys- 
tem would be a 'nominal' band insulator even in the ab- 
sence of any antiferromagnetic order. In 2D the ground 
state is known [|l8|,^ to have antiferromagnetic order, 
which in principle could turn even the unhybridized con- 
duction electron system into an insulator, due to the re- 
duction of the Brillouin zone by a factor of two. However, 
our data do not show any indication of antiferromagnetic 
ordering in the c-like single-particle spectrum or the c- 
like momentum distribution, so that the 'effective' SDW- 
potential felt by the c-electrons due to the ordering of the 
/-electrons must be extremely weak. Moreover, the band 
structure obtained from simple SDW mean field calcula- 
tion is qualitatively different from the numerical results, 
and in particular would lack the dispersionless low en- 
ergy band with /-character, which are typical of the band 
structures in both ID and 2D. By analogy with the ID 
system we thus conclude that antifcrromagnctism is not 
essential for the insulating nature of the ground state 
(as can be seen also from the fact that the strong cou- 
pling model remains an insulator in the spin-liquid phase 
pO|), but that it is the 'merging' of the c- and /-electron 
systems which drives the metal-insulator transition with 
decreasing temperature. 

Increasing the temperature beyond Tcoh the /-electrons 
drop out of the Fermi surface volume, which results in the 
metal-insulator transition. The /-electrons do no longer 
participate in the low energy physics, the first ionization 
states are purely c-like and the /-spin excitation becomes 
localized. In the single-particle spectrum we can iden- 
tify dispersionless /-like sidebands, well separated from 
the chemical potential. This suggests that local singlets 
between / and c-electrons still exist, but no phase co- 
herence between the singlets is established. Finally, at 
the highest temperatures the /-electrons disappear com- 
pletely from the low energy single-particle spectrum, and 
all that remains are two /-like Hubbard bands and the 
standard free-electron band for the c-electrons. We as- 
sociate the temperature where this complete decoupling 
occurs as the analogue of the Kondo temperature. 
A surprising feature of the results is the close similarity 
with those obtained for the ID system In particular, 
the characteristic temperatures are more or less identical 



(to the accuracy to which these 'crossover' temperatures 
can be assigned) in ID and 2D. Together with the very 
different nature of the single-particle DOS around the 
band center in the ID and 2D tight-binding bands this 
seems to indicate that in the lattice case there is no sig- 
nificant dependence of the characteristic temperatures on 
the density of states at the Fermi energy. 
At low temperatures our data in 2D are consistent with 
long range antiferromagnetic order, in agreement with 
previous work by Vekic et al. and by Assaad . How- 
ever, our data do not show any indication of 'antifer- 
romagnetic symmetry' in the single particle spectra, in 
particular there are no distinguishable antiferromagnetic 
umklap bands in the spectra. This again indicates that 
an SDW-like treatment of the antiferromagnetic state is 
not adequate. 
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